The study of leading twist light cone wave function of rjc meson. 
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This paper is devoted to the study of leading twist hght cone wave function of r/c meson. The 
moments of this wave function have been calculated within three approaches: potential models, 
nonrelativistic QCD and QCD sum rules. Using the results obtained within these approaches the 
model for the light cone wave function of leading twist has been proposed. Being scale dependent 
light cone wave function has very interesting properties at scales ^ > mc'. improvement of the 
accuracy of the model, appearance of relativistic tail and violation of nonrelativistic QCD velocity 
scaling rules. The last two properties are the properties of real leading twist light cone wave function 
of rjc meson. 

PACS numbers: 12.38.-t, 12.38.Bx, 13. 66. Be, 13.25. Gv 

I. INTRODUCTION 



Commonly exclusive charmonium production at high energies is studied within nonrelativistic QCD (NRQCD) fl*]. 
, In the framework of this approach charmonium is considered as a bound state of quark-antiquark pair moving with 
• small relative velocity v 1. Due to the presence of small parameter v the amplitude of charmonium production can 
I be built as an expansion in relative velocity v. 

. Thus in the framework NRQCD the amplitude of any process is a series in relative velocity v. Usually, in the most 
of applications of NRQCD, the consideration is restricted by the leading order approximation in relative velocity. 
However, this approximation has two problems which make it unreliable. The first problem is connected with rather 
large value of relative velocity for charmonium v'^ ~ 0.3, v ~ 0.5. For this value of one can expect large contribution 
, from relativistic corrections in any process. So in any process resummation of relativistic corrections should be done or 
' one should prove that resummation of all terms is not crucial. The second problem is connected with QCD radiative 
T— I I corrections. The point is that due to the presence of large energy scale Q there appear large logarithmic terms 
■ (tts log Q/toc)", Q 3> rric which can be even more important than relativistic corrections at sufficiently large energy 
( Q ^ 10 GeV). So this terms should also be resummed. In principle, it is possible to resum large logarithms in the 
NRQCD factorization framework 0, [3], however such resummation is done rarely. 

The illustration of all mentioned facts is the process of double charmonium production in e^e~ annihilation at 
Oh, B-factories, where leading order NRQCD predictions [A,, Jj] are approximately by an order of magnitude less than 
Q_^' experimental results 0,[8[- The calculation of QCD radiative corrections [§) diminished this disagreement but did not 
' remove it. Probably the agreement with the experiments can be achieved if, in addition to QCD radiative corrections, 
relativistic corrections will be resummed fTo| . 
^ ' In addition to NRQCD, hard exclusive processes can be studied in the framework of light cone expansion formalism 
• ^ , [Til [T^ where both problems mentioned above can be solved. Within light cone expansion formalism the amplitude is 
^ ' built as an expansion over inverse powers of characteristic energy of the process. Usually this approach is successfully 
^ , applied for excusive production of light mesons T?| . However recently the application of light cone expansion formalism 
- - ' to double charmonium production ^isl . 0, [l5l . ,16.] allowed one to achieve good agreement with the experiments. 

In the framework of light cone formalism the amplitude of some meson production in any hard process can be 
written as a convolution of the hard part of the process, which can be calculated using perturbative QCD, and process 
independent light cone wave function (LCWF) of this meson that parameterizes nonperturbative effects. From this one 
can conclude that charmonium LCWFs are key ingredient of any hard exclusive process with charmonium production. 
This paper is devoted to the study of leading twist LCWF of rjc meson. 

The paper is organized as follows. In the next section all definitions needed in our calculation will be given. In 
Section HI the moments of LCWF will be calculated in the framework of BuchmuUer-Tye and Cornell potential 
models. Section VI is devoted to the calculation of the moments within NRQCD. QCD sum rules will be applied to 
the calculation of the moments in Section V. Using the results obtained in Sections III-V the model for LCWF will 
be built in Section VI. In Section VII this model will be compared with some other models proposed in literature. In 
the last section we summarize the results of our paper. 
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II. DEFINITIONS. 

The leading twist light cone wave function (LCWF) of r/c meson can be defined as follows 

{0\Q{z)j^J5[z,-z]Q{~z)\^,)^=zUp^ I ^d^e'^P'^^<l>{^,fi), (1) 

where the following designations are used: xi , X2 are the parts of momentum of the whole meson carried by quark 
and antiquark correspondingly, ^ = xi — X2, p is a momentum of rjc meson, fi is an energy scale. The factor [z, —2], 
that makes matrix element ([T]) gauge invariant, is defined as 

[z,-z] = PeM^9 J dxt'A^ix)]. (2) 

The LCWF /i) is normalized as 

= (3) 



-1 

With this normalization condition the constant /r,^ is defined as 

(0|Q(0)7a75Q(0)|r/c> (4) 
LCWF (/)(a;,/i) can be expanded in Gegenbauer polynomials Cn^^iO follows 
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1+ J2 ^n{p^)Cl''{i) 
n=2,4.. 



(5) 



At leading logarithmic accuracy the coefficients a„ are renormalized multiplicatively 



""(m) = ( "''/^l ) an{^J-o), (6) 



where 



4/ 2 . 2 



1 \ 2 
4^- , 6o = ll--nfl. (7) 



3V (n + l)(n + 2) f^jj' 3 

It should be noted here that conformal expansion ([5]) is a solution of Bethe-Salpeter equation with one gluon exchange 
kernel 

From equations ©-([T]) it is not difficult to see that at infinitely large energy scale — > 00 LCWF (f>{^, /i) tends to 
the asymptotic form <j)as{£.) — 3/4(1 — ^^). But at energy scales accessible at current experiments the LCWF 0(f , fi) 
is far from it's asymptotic form. The main goal of this paper is to calculate the LCWF /i) of r/c meson. One 
way to do this is to calculate the coefficients of conformal expansion ^ an- Having this coefficients at some energy 
scale ^0 one can find the LCWF (/)(^, fi) at any scale /i. Alternative method to parameterize LCWF is to evaluate it's 
moments at some scale 

(Dm- j\d£,C ^^). (8) 

In our paper this method will be used. It is worth noting that since 7]c mesons have positive charge parity the LCWF 
(j){£,,fj,) is ^-even. Thus all odd moments (^2fc+i^ equal zero and one needs to calculate only even moments. 

Below we will need the formula that connects moment (^") with the matrix element {0\Qj^j5{iz°'Da-)"Q\P{p))- 
To obtain it we expand both sides of equation ([1]) 



n ' n ' 

and get 

(0|Q7.75(*^"i3,)"g|77c) = ^Up.izpnO- (10) 

Here 

D = D-D, D = 1) -igB^iX" /2). (11) 
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III. THE MOMENTS IN THE FRAMEWORK OF POTENTIAL MODELS. 

In potential models charmonium mesons are considered as a quark- antiquark system bounded by some static 
potential. These models allow one to understand many properties of chamonium mesons. For instance, the spectrum 
of charmonium family can be very well reproduced in the framework of potential models [iTj . Due to this success one 
can hope that potential models can be applied to the calculation of charmonium equal time wave functions. 

Having equal time wave function of 77^ meson in momentum space one can apply Brodsky- Huang-Lepage 

(BHL) [l^ procedure and get the LCWF of leading twist fi) using the following rule: 

Htf^)- / cPk^Mx,^±), (12) 



where 'ipc{x,'k±) can be obtained from ipCk) after the substitution [T^l 

k^-k^, h^{x,~x,)^, M',=^^l±^. (13) 

2 X1X2 

Here Mc is a quark mass in potential model. In our paper equal time wave function V)(k) will be calculated in 
the framework of the potential models with BuchmuUer-Tye pjj] and Cornell potentials [2Q|- The parameters of 
BuchmuUer-Tye potential model were taken from paper [l^. For Cornell potential V{r) ~ —k/r + r/a? the calculation 
was carried our with the following set_of parameters: k = 0.61, a = 2.38 GeV~^, Mc — 1.84 GeV }21l] . 

It is worth noting that in paper ^22 | the relations between the light cone wave functions and equal time wave 
functions of charmonium mesons in the rest frame were derived. The procedure proposed in paper [l^l is similar 
to BHL with the difference: in formula one must make the substitution cP k±^yW + rn^/ {Am cX 1X2). 

But this substitution was derived at leading order approximation in relative velocity of quark-antiquark motion 
inside the charmonium. At this approximation ~ 0{v'^), AxiX2 ~ 1 -|- 0{v^) and the substitution amounts to 
(Pk^ (Pk±{l + 0{v^)). Thus at leading order approximation applied in [2^ these two approaches coincide. 

The potential models with BuchmuUer-Tye and Cornell potentials treat charmonium as a nonrelativistic bound 
state, so these models can not be applied in the region where the motion is relativistic. From this one can conclude 
that it is not possible to apply expression p2)) at sufficiently large scale /i ^ Mc- At the same time conformal 
expansion ([5]) can be applied at scales /i^ ^ ~ M^v^ where q and v are characteristic momentum and velocity of 
cc system. Thus to find the values of moments (^") one should apply formula (fT2|l at scale not too low and not too 
large. Our calculation will be done at scale ^ = 1.5 GeV ~ Mc where both conformal expansion ([5]) and formula 
can be applied. To estimate the scale dependence of our results in addition to scale ^ = 1.5 GeV we have also done 
the calculation at scale fj, = 2 GeV. The moments at scale = 2.0 GeV differ by few percents from the moments 
calculated at scale = 1.5 GeV so the scale dependence is very weak. 

The results of our calculation are presented in Table I. In second and third columns the moments calculated in 
the framework of BuchmuUer-Tye and Cornell models are presented. It is seen that there is good agreement between 
these two models. 

It should be noted here that the larger the power of the moment the larger contribution form the end point regions 
(x ~ and a; ~ 1) it gets. From formulas (|13p one sees that the motion of quark-antiquark pair in these regions is 
relativistic and cannot be considered reliably in the framework of potential models. Thus it is not possible to calculate 
higher moments within the potential models. Due to this fact we have restricted our calculation by few first moments. 

IV. THE MOMENTS IN THE FRAMEWORK OF NRQCD. 

The moments of LCWF ([T]) can be calculated in the framework of NRQCD. To do this let us consider matrix 

element {0\Q"f^"f^{iD^-^)(iD^^)Q\ric) in the meson's center mass frame. In this section we will work at leading order 
approximation in relative velocity of quark-antiquark pair v. According to the velocity scaling rules gauge-covariant 

time derivative Dt is suppressed by one power of v in comparison with gauge-covariant spatial derivative Di. Thus 
if we set fii or ^2 to zero than at leading order approximation this matrix element is zero. Moreover in the meson's 
center mass frame only zero component of the current v — differs from zero. Regarding these two properties it 
causes no difBculties to built the matrix element: 

{0\Q-f^-f5{iD^^){iD^^)Q\ric) = ifp^T^^^^ + higher order in v terms, (14) 



4 



where tensor T^^^jj is defined as 



Vc 



The constant / can be calculated contracting the indices /ii and ^2 and regarding the fact that time derivative is 
suppressed 

/=^(0|Q7o75(«5fQ|ry,). (16) 

The constant / can also be related to the moment (^^) of the LCWF. To obtain this relation let us consider matrix 
element (1141) in the infinite momentum frame {p — > 00). Obviously only tensor structure p^iPfj.2 gives contribution to 
the LCWF of leading twist (H]) . At the same time at leading twist accuracy matrix element (ITil) can be found from 
equation pO|) . Thus we get 

(17) 

Combining equations (fTB)) . (fTZ)) and the definition of the constant /^^ ^ we obtain 



1 (0|Q7o75(»D)^Q |^c 
SMI (0|Q7o75Qhc 



) = ¥7?2 /ni/n.^ ^1^ N — • (18) 



At leading order approximation M,,^ = 2Mc 



(0|g7o75(*D)2g|r;,) ^ 4(0|x+(*D)2^|r;,), (19) 
(0|g7o75g|r/c> - (Olx+V-hc), (20) 

where V' and are Pauli spinor fields that annihilate a quark and an antiquark respectively. The factor 4 in equation 
()19p appears since the distance between quark and antiquark in left hand side of this equation is 2z and z in the right 
hand side. With these relations equation can be rewritten as 

^ 3M2 (0|x+^|77c) 3 ■ ^ ^ 

Analogously it is not difficult to calculate similar expressions for higher moments. We will write the expressions for 
{^*) and (^^) without derivation 

4 ^ 1 (0|x+(^D)V|r/c) ^ , . 

^ 5Af4 (0|x+^|r7e> 5 ' ^ ' 

(0|x+(»D)Vl^c) ^ {j^ , . 

^ 7Af6 (0|x+V^|?7c) 7 ■ ^ ^ 

It should be noted here that equations (|2ip - (|23p were derived within dimensional regularization at leading order 
approximation in as- Moreover, in equations (pT|) - (p3| the moments are defined at scale ~ Mc since NRQCD 

matrix elements (0|x^V'l'7c), (0|x^(*D)^V'|'7c)7 ■•• that will be used in our calculations [2^ were derived at this scale. 

In the framework of NRQCD the matrix elements (0|x^(iD)"i/'|77c) were calculated in paper [2^ 

where Eb is a bound state energy of rjc meson, m is the mass of quark in paper [2^. The main idea of paper [2^ is 
that at leading order approximation in relative velocity NRQCD results can be reproduced by potential model with 
heavy quark potential measured in lattice. The results of the calculation of the moments are presented in the fourth 
column of Table I. In the evaluation of the moments we have used the following values of parameters: = M^^/2, 
the mass ra and bound state energy Eb are taken from paper (23| . The central values of the moments and the 
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BuchmuUer-Tye 
model [19J 


Cornell 
model [201 


NRQCD 
[23] 


QCD 
sum rules 


n = 2 


0.086 


0.084 


0.075 ±0.011 


0.070 ± 0.007 


n = 4 


0.020 


0.019 


0.010 ± 0.003 


0.012 ±0.002 


n = 6 


0.0066 


0.0066 


0.0017 ± 0.0007 


0.0032 ± 0.0009 



TABLE I: The moments of LCWF obtained within different approaches. In the second and third columns the moments 
calculated in the framework of BuchmuUer-Tye and Cornell potential models are presented. In the fourth column NRQCD 
predictions for the moments are presented. In last column the results obtained within QCD sum rules are shown. 



errors have been calculated according to the approach proposed in the same paper. It should be noted that the error 
presented in Table I is assigned to the uncertainty of potential model parameters used in the calculation. In addition 
to this error there is an uncertainty due to higher order v corrections. For the second moment one can expect that 
the error is about 30%. For higher moments this error is larger. 

It is seen from Table I that NRQCD predictions for the second and the fourth moments are in good agreement 
with potential model and there is disagreement for the moment (^^) between these two approaches. The cause of 
this disagreement is the fact noted above: due to the large contribution of relativistic motion of quark-antiquark pair 
inside quarkonium it is not possible to apply both approaches for higher moments. We believe that both approaches 
can be used for the estimation of the values of the second and the fourth moments. The predictions for the sixth and 
higher moments obtained within both approaches become unreliable. 

The aim of this paper is the study of leading twist LCWF of rjc meson. But at leading order approximation in 
relative velocity v there is no difference between rjc meson and J/^* mesons. So the results for the moments obtained 
within potential models and NRQCD are valid for leading twist LCWFs of J/^ meson. Moreover, at leading order 
approximation in relative velocity the moments of leading twist LCWF of rjc meson equal to the moments of higher twist 

LCWFs of rjc and J/^* mesons since these moments can be expresses through the matrix elements: (0|x+(iD)^-!/)|77c), 
(0|x+(*D)Vl^c), ... 



V. THE MOMENTS IN THE FRAMEWORK OF QCD SUM RULES. 

In the previous sections two approaches to the calculation of the moments of leading twist LCWF were considered. 
The main disadvantage of both approaches is that quarkonium is considered as a nonrelativistic bound state of quark- 
antiquark pair. Due to this the errors of the calculation for both approaches are rather large. Another approach to 
the study of the moments of LCWF based on QCD sum rules [U, |25[ was developed by Chernyak and Zhitnitsky in 
paper [26i]. This approach does not assume that the quarks in quarkonium are nonrelativistic. In this section QCD 
sum rules will be applied to the calculation of the moments of LCWF. 

To calculate the moments of LCWF in the framework of QCD sum rules let us consider two-point correlator: 

Ujz,q) - z J d^xe''^^Q\TJo{x)JnmO) = {zqr+^Uniq^), (25) 

Mx) = Q{x)z^5Qix), J„(0) = Q(0)l75(iz'^S„)"Q(0), = 0. 

It should be noted that explicit dependence of QCD correlator Iln{q^) on n will be dropped in subsequent equations. 

The calculation of QCD expression for two-point correlator TlQcuiQ^) is done through the use of operator product 
expansion (OPE) for the T-ordered product of currents. As a result of OPE one obtains besides usual perturbative 
contribution TlpcitiQ^) a-lso nonperturbative power corrections nnpoit('?^), given by QCD vacuum condensates. The 
diagram that gives the leading order perturbative contribution is shown in Fig. la. The leading order nonperturbative 
power corrections are given by gluon vacuum condensate (G^i/G^''). The diagrams that contribute to correlator ([^5]) 
are shown Fig. Ib-d. The connection to hadrons in the framework of QCD sum rules is obtained by matching the 
resulting QCD expression for current correlator with spectral representation, following the structure of dispersion 
relation at g'^ < 0: 

nphys(9 ) = / ds — + subtractions, (26) 

J4m^ S q 
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abed 

Fig.l: The diagrams that contribute to the perturbative (a) and nonperturbative (b,c,d) parts of QCD sum rules. 

Assuming that the dispersion relation (j26p is well convergent, the physical spectral functions are generally saturated 
by the lowest lying hadronic states plus a continuum starting at some thresholds s„: 



(27) 



where 



(28) 



Continuum contribution p'^o"* for s > s„ can be approximated by spectral density of perturbative contribution 
~ Imnpert(Q^)- To get QCD sum rules for the moments of LCWF one should apply operation 



(29) 



to both QCD and physical expressions for correlator ((25|) and then equate them. Thus we get QCD sum rules for the 
moments: 



1 



(ml + g2)m+l ^ J^^, (s + Q2)m+1 



ds 



imn 



port V 



(30) 



where Imnport(s), ^^J^lrtiQ^) can be written as 



Imnpert(s) = A^;"+l(^ - 

OTT n + 1 n + 6 



Ami 



(31) 



n 



(™) (C)2 
npcrtV^ 



n(")(Q2) 
n('")(Q2) 



n("'(Q2) 



(l_^2)m+2 



2(m2 + 3m + 2) / ^^^(1 + 3^') 



(4m2 + Q2(i_^2)) 

(1_^2)™+1 



m+2 ■ 



(4to2+Q2(1_^2)) 



m+3 ' 



'-{n^ -n)(m+l) df C~ X9, 

384^^ V-l ' (4^2+Q2(l_^2))"+2' 



(32) 



here = —q^- In the original paper [1^ the method QCD sum rules was applied at ~ 0. However as it was 
shown in paper p7j there is large contribution of higher dimensional operators at — q which grows rapidly with 
TO. To suppress this contribution in our paper sum rules (j30p will be applied at = 4to2. 

In the numerical analysis of QCD sum rules the values of parameters and (a^G^/Tr) will be taken from paper: 

M 



mc 



1.24 ± 0.02 GeV, (— G^) = 0.012 ± 30% GeV*. 



(33) 



First sum rules ([3D]) will be applied to the calculation of the constant /^^ . It is not difficult to express the constant /^^ 
from equation (j30p at n = as a function of m. For too small values of m (to < toi) there is large contributions from 
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higher resonances and continuum which spoil sum rules (j30p . Although for m ^ mi these contributions are strongly 
suppressed, it is not possible to apply sum rules for too large m (m > 7712) since the contribution arising from higher 
dimensional vacuum condesates rapidly grows with m what invalidates our approximation. If mi < m2 there is some 
region of applicability of sum rules pop [mi, m2] where the resonance and the higher dimensional vacuum condensates 
contributions arc not too large. Within this region /^^ as a function of m varies slowly and one can determine the 
value of this constant. The value of the continuum threshold s„ must be taken so that to appear stability region 
[mi, m2]. Our calculation shows that for central values of the parameters p3p there is stability in the region m > 4 if 
one takes s„ > 3.6^ GeV^. In paper [28], where the mass of rjc meson was studied within QCD sum rules, the authors 
took Sn = 00. In our paper we will use the same value of the threshold. Applying approach described above we get 

= 0.120 ± 0.002 ± 0.005 ± 0.016 GeV^ (34) 

The region of stability in the determination of the constant /^^ begins at mi = 6 and ends at m2 = 10. The first 
error in (j34p corresponds the variation of /^^ within this region. The second and the third errors in (I34p correspond 
to the variation of the gluon condensate {agC^) and the mass rric within ranges p3)) . From the results ((34)) one sees 
that the main error in determination of the constant /^^ results from the variation of the parameter rric. This fact 
represents well known property: high sensitivity of QCD sum rules to the mass parameter rric. 

Next let us consider the second moment of LCWF (<^^) in the framework of QCD sum rules. One way to find the 
value of (^^) is to determine the value of from sum rules ((30)) at n = 2 and then extract (^^). However, as 

it was noted above, this approach suffers from high sensitivity of right side of equation (|30p to the variation of the 
parameter mj,. Moreover, the quantity /^(C^) includes not only the error in determination of (^^), but also the error 
in To remove this disadvantages in our calculation we will consider the ratio of sum rules at n = 2 and n — 0: 
Irii^"^) I fri- The moments (^^), (^®) will be considered analogously. Applying standard procedure one gets the results: 

(^2) ^ 0.070 ± 0.002 ±0.007 ±0.003, (35) 
{S,^) = 0.012 ±0.001 ± 0.002 ±0.001, 
{S,^) = 0.0032 ±0.0002 ±0.0009 ±0.0003. 

In the calculation of moments the stability was achieved within the interval m e [9,13]. The first error in ()35p 
corresponds the variation within the region of stability. The second and the third errors in (I35p correspond to the 
variation of the gluon condensate (asG^) and the mass nic within ranges ()33p . It is seen that, as one expected, the 
sensitivity of the ratio /,^(^") //^ to the variation of rric is rather low. The main source of uncertainty is the variation 
of gluon condensate (asC^). In the fourth column of Table I the results of our calculation are presented. The errors 
in Table I correspond to the main source of uncertainty — the variation of gluon condensate (asC^). 

It is not difficult to show that in the calculation of QCD expression for the correlator ((25)) characteristic virtuality 
of quark is ~ (4m^ + Q^)/m ~ m^. So the values of the moments ()35p are defined at scale ~ m^. It should be noted 
here that, in all approaches applied to the calculation of the moments, the quantities (^") are defined at different 
scales. But the scale dependence appears at NLO as correction and since all these scales are of order of ~ nic, the 
variation of the scale near ~ rric will not change our results considerably. For this reason we believe that it is possible 
to compare the results obtained within different approaches. 

From Table I it is seen that the larger the number of the moment n the larger the uncertainty due to the variation 
of vacuum gluon condensate. This property is a consequence of the fact that the role of power corrections in the 
sum rules pop grows with n. From this one can conclude that there is considerable nonperturbative contribution to 
the moments (^") with large n what means that nonperturbative effects are very important in relativistic motion of 
quark-antiquark pair inside the meson. The second important contribution to QCD sum rules ()30p at large n is QCD 
radiative corrections to perturbative part npcrt(Q^)- Unfortunately today one does not know the expression for these 
corrections and for this reason we cannot include them to sum sules (I30p . We can only say that these corrections 
grow with n and probably the size of radiative corrections to the ratio /^(C")//^ is not too big for not too large n. 
Thus one can expect that QCD radiative corrections will not change dramatically our results for the moments n = 2 
and n = 4. But the radiative corrections to (f^) may be important. 

It is interesting to estimate relative velocity of quark-antiquark pair inside rjc meson using the results obtained 
within QCD sum rules. In the section devoted to the calculation of the moments in the framework of NRQCD we 
have derived formula ()2ip . At leading order approximation this formula allows one to connect the square of relative 
velocity and the second moment of LCWF. Applying this formula one gets 

{v^)r,, = 3(^^),,^ = 0.21 ± 0.02. (36) 

This value is in good agreement with potential models estimations of (w^) [l^, [2^. The error in ((55)) corresponds to 
the main source of uncertainty of our results — the variation of gluon condensate (asG^). In addition to this error 
there is an uncertainty due to higher order v corrections ~ 30%. 
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Fig.2: The LCWF (|38j at scales /io = 1.2 GeV, /ii = 10 GeV, ^2 = 100 GeV,^3 = 00. 

It should be noted here that in our calculation we have used the value of the threshold s„ = cx). But if one takes 
the value of the threshold energy within a reasonable interval this will lead to the shift of the stability region to lower 
m and the results of our calculation will be changed not greater than by few percents for (^^), (^^) and 10% for (^^). 

We believe that the values of the moments ([55]) calculated within QCD sum rules are more reliable than the ones 
obtained in the framework of potential models and NRQCD. For this reason we consider results (p5|) as the main 
result of our paper and below these values will be used. 



VI. THE MODEL FOR THE LCWF OF ry^ MESON. 

To build the model of the leading twist LCWF of ijc meson let us write Borel version [l^, [HI of sum rules ([50]) 
without continuum contribution p'^o"* and power corrections Ilnport 



fliC)e-<J-'^^ .er,(l-aexp -^^-^ . (37) 



Within this approximation the LCWF can be written in the form 

0(e, fi^m,) = c(/3)(l - C')exp (" ^ (38) 

where c(/3) is a normalization constant and /3 is a constant that will be calculated below. Although sum rules (|37l) is 
much simpler that sum rules pop . to reproduce the results obtained in the previous section with good accuracy it is 
sufficient to use LCWF in the form of equation (|38p. To fix the constant f3 the moment (^^) will be taken. Moreover 
we will suppose that the results (P5|) obtained in the previous section are defined at scale ^0 = 1.2 GeV ^ rric- Thus 
we get /? = 3.8 ± 0.7. The constant c(/3) can be determined from normalization condition ([3]). It is not difficult to see 
that 1/ (3 ^ as one can expect. The moments of this wave function are 

(^^) = 0.070 ± 0.007, (39) 
(^4) = 0.012 ± 0.002, 
(^^) = 0.0030 ± 0.0009. 

At central value /3 ~ 3.8 the constant c{(3) ~ 62. Below the LCWF ([38| at central value will be used. 

LCWF ([55]) is defined at scale fi = fig. It is not difficult to calculate this function at any scale fi > fiQ using 
conformal expansion The LCWFs at scales /Xq — 1.2 GeV, /^i = 10 GeV, 112 = 100 GeV, 113 = 00 are shown in Fig. 
2. The moments of LCWFs at scales no = 1.2 GeV, ^1 — 10 GeV, /Z2 = 100 GeV, ^3 = 00 are presented in second, 
third, fourth and fifth columns of Table II. 

The model p8p has some interesting properties. From conformal expansion ([5]) one can derive the expressions that 
determine the evolution of the moments: 

(^'>M = ^+a2(A*)§, (40) 

- |+a2(A^)4+a4(A^)^, 
/^6\ 1 / ^12 / N 120 , , 64 
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Similar relations can be found for any moment. Further let us consider the expression for the second moment (^^). 
In our paper we have found the value (^^) with some error at scale /i = /ig- This means that the value of the 
coefficient a2(/i = fJ-o) was found with some error. The coefficient 02 decreases as scale increases. So the error in 
02 and consequently in (^^) decreases as scale increases. At infinitely large scale there is no error in (^^) at all. 
Our calculations show that the error 10% in (^^) at scale jj, — hq decreases to 4% at scale /i = 10 GeV. Applying 
relations (|40|) it is not difficult to show that similar improvement of the accuracy happens for higher moments. The 
improvement of the accuracy allows us to expect that the model ([55)1 at larger scales will be rather good even if QCD 
radiative corrections to the results (|35p are large. 

From Fig. 2 one sees that LCWF at scale fi — fiQ practically vanishes in the regions 0.75 < |^| < 1. In this region 
the motion of quark-antiquark pair is relativistic and vanishing of LCWF in this region means that at scale /x = /xg 
charmonium can be considered as a nonrelativistic bound state of quark-antiquark pair with characteristic velocity 
1//3 ~ 0.3. Further let us regard the function — /ig) as a conformal expansion ([5]). To get considerable 
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suppression of the LCWF in the region 0.75 < |^| < 1 one should require fine tuning of the coefficients of conformal 
expansion a„(/i = ^g). The evolution of the constants a„( especially with large n) near /i = /ig is rather rapid ( see 
formulas ^ and ([7])) and if there is fine tuning of the constants at scale /i = /xg this fine tuning will be rapidly broken 
at larger scales. This property is well seen in Table II and Fig. 2. From Fig. 2 it is seen that there is relativistic tail 
in the region 0.75 < |^| < 1 for scales = 10, 100 GeV which is absent at scale fi = fiQ. Evidently this tail cannot 
be regarded in the framework of NRQCD. This means that, strictly speaking, at some scale charmonium can not 
be considered as nonrelativistic particle and the application of NRQCD to the production of charmonium at large 
scales may lead to large error. Although in our arguments we have used the model of LCWF ([55]) it is not difficult to 
understand that the main conclusion is model independent. 

According to the velocity scaling rule [l| the moments (^") of LCWF depend on relative velocity as ~ u". It is not 
difficult to show that the moments of LCWF satisfy these rules. Now let us consider the expressions that allows 
one to connect the coefficients of conformal expansion a„ with the moments (^"). These expressions for the moments 
(C''): (^^) are given by formulas (|40|1 . It causes no difficulties to find similar expressions for any moment. From 
expressions (|40l) one sees that to get velocity scaling rules: (^") ~ v" at some scale one should require fine tuning of 
the coefficients a„ at this scale. But, as was already noted above, if there is fine tuning of the coefficients at some 
scale this fine tuning will be broken at larger scales. From this one can conclude that velocity scaling rules are broken 
at large scales. 

Consider the moments of LCWF ([55)1 at infinite scale. It is not difficult to find that 

(r)M = oo = , ^^y (41) 

[n + l){n + 3) 

From last equation one can find that (^") does not scale as as velocity scale rules [l| require. Thus scaling rules 
obtained in paper [1] are broken for asymptotic function. Actually one does not need to set the scale fi to infinity to 
break these rules. For any scale /i > ^g there is a number ng for which the moments (^"'), n > uq violate velocity 
scaling rules. This property is a consequence of the following fact: beginning from some n = ng the contribution of 
the relativistic tail of LCWF, that appears at scales /i > /xg, to the moments becomes considerable. 

The amplitude T of any hard process with charmonium meson production can be written as a convolution of LCWF 
<I>(^) with hard kernel i?(^) of the process. If one expands this kernel over ^ and substitute this expansion to the 
amplitude T one gets the results: 

If one takes the scale /i ~ /ig in formula (|^^ than moments (^") scale according to the velocity scaling rules ~ w" 
and one gets usual NRQCD expansion of the amplitude. However due to the presence of the scale of the hard 
process fih ^ Mo there appears large logarithms log/i^Z/ig which spoil NRQCD expansion (|42|) . To remove this large 
logarithms one should take /i ~ /i/j. But at large scales velocity scaling rules are broken and application of NRQCD 
is questionable. 



VII. COMPARISON WITH OTHER MODELS. 



Asymptotic wave function. 

We have already noted that if in conformal expansion ([S]) the scale /i tends to infinity LCWF tends to the function 

0a. = ^(1- a, (43) 
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<^(C,Mo = 1.2 GeV) 


<j>{^,^Il = 10 GeV) 


(j>{^, fj,2 = 100 GeV) 


= OO) 


BC [13] 


BKL 122} 


n = 2 


0.070 


0.12 


0.14 


0.20 


0.13 


0.019 


n = 4 


0.012 


0.040 


0.052 


0.086 


0.040 


0.0083 


n = 6 


0.0032 


0.019 


0.026 


0.048 


0.018 


0.0026 



TABLE II: The moments of LCWF ((351) proposed in this paper at scales fio = 1.2 GeV,/ii = 10 GeV /i2 = 100 GeV, ^3 = 00 
are presented in second, third, fourth and fifth columns. The moments of LCWF proposed in paper [1^ are shown in the sixth 
column. The moments of LCWF proposed in paper (2^] are shown in the sixth column. 



which is called asymptotic. This wave function does not depend on the mesons with the same quantum numbers. For 
instance, there is one LCWF defined by equation (|43p for r/c and tt mesons. From this point of view the distinction of 
quarkonium LCWF from it's asymptotic form can be considered as a measure of relativism of the quarkonium. For 
particularly nonrelativistic object the relation (^^) ^ {£,'^)fj.=oc must be satisfied. 

The moments of asymptotic LCWF are presented in the fifth column of Table II. One sees that although at scale 
M = Mo '7c can be considered as nonrelativistic object with v'^ ~ 0.3, the relation for the particularly nonrelativistic 
object (^^)^=oo ^ {£,'^)fj.=m^ is badly satisfied. At scale /i = 10 GeV this relation is not satisfied at all and at this 
scale ?7c can not be considered as a nonrelativistic object. This fact has already been noted in the previous section. 

The LCWF proposed by Bondar and Chernyak in paper |l3l |. 

To resolve contradiction between NRQCD prediction Q and experimental measurements of the cross section of the 
process e+e" — > J/'^rjc at energy y^s = 10.6 GeV 0,3 Bondar and Chernyak proposed the following model of LCWF 

^Bc(.0-civ^)<Pasi0 L, . ""'^I 2nt| ' (44) 

[1 - 4X1X2(1 - v^)\ J 

where = 0.3 is relative velocity of quark-antiquark pair in ijc meson, c(v'^) is the coefficient which is fixed by the 
wave function normalization J dx(j)Bc{C) — 1> 0as(C) is the asymptotic LCWF (|43| . It should be noted that LCWF 
breaks velocity scaling rules. 

The moments of this wave function are presented in the sixth column of Table II. According to paper [l^ LCWF 
is defined at scale ^ ^ hq- However from Table II it can be seen that the second moment of this function is 
approximately two times larger than the second moments calculated in this paper. For higher moments the difference 
is greater. The moments of this wave function are in better agreement with the moments of — 10 GcV). Our 

calculation shows that within the region ^ e (—0.8,0.8) the function = 10 GeV) differs from function by 

less than 15%. In the end point region |^| ~ 1 these functions considerably differ from each other but this region does 
not give main contribution to the cross section of the process e~^e~ J/'^rjc- The fact that the scale /i = 10 GeV is 
an energy at which the cross section of the process e"'"e~ — > J/'^rjc was measured 0, [1] allows us to understand why 
application of LCWF led to a good agreement with the experiments. 

Effective LCWF proposed by Bodwin, Kang and Lee (BKL) in paper |22| . 

To reconcile light cone and NRQCD approaches to calculating e+e^ J/'if + ric the authors of paper [22*1 proposed 
LCWF of rjc meson. In our paper we are not going to discuss the approach applied in paper ^22^ . Only main results 
of this paper will be considered. Thus LCWF obtained in this paper is 




dp \/ + ml I dr ip{r) sin(pr), (45) 
where 

diO=ml^^, (46) 

TOc is the pole mass of c quark, ipi''^) is a solution of Schrodinger equation with the potential —k/r + rjo?. In our 
calculation potential model parameters and pole mass are taken from paper [2^ . LCWF (|45)) will not be discussed 
here since LCWF (j45|) and (|44p are similar in shape and have approximately equal moments. It should be noted here 
only that LCWF is wider than model of LCWF ([55)1 proposed in this paper. 

Expression (I45p contains high relative momentum-p tail which cannot be regarded within nonrelativistic potential 
model. To resolve this problem the authors of paper j22| introduced effective LCWF that removed this high momentum 
tail from LCWF (HI). 
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The resulting effective LCWF is 

4>BKL{0=M0~^m, (47) 

where 



Sir rrir 



e 2 "i-yr^ 



(48) 



In the last equation the following designations were used: m and is the mass of c quark and coupling constant in 
potential model [l^l and [/(C)]+4 means 



£^ d^f ioumo = duiomo - Hio)]. (49) 

As it was noted above effective function (|47p is not true LCWF. The subtraction A(/)(^) is a formal device to remove 
the high-momentum contributions that are not calculated reliably in a potential model and to avoid double count 
the NRQCD order-as corrections. However, we believe that one can compare the LCWF (|38p and function ([TT]) 
for the following reasons. First, as it was noted above at scale /i ~ /io there is no high-momentum contribution in 
model ([55)1 . Moreover, the calculation of the moments of LCWF was done within QCD sum rules at leading order 
approximation in a^, so there is no double counting in the framework of light cone formalism^. Second, the main 
goal of this comparison is to understand what results one can expect if instead of LCWF ([38|) function (|47|) will be 
applied. 

The moments of effective LCWF (l47l) are presented in the seventh column of Table II. From this table it is seen 
that the moments of model (|47p are much smaller not only than the predictions obtained within potential models and 
QCD sum rules but also than NRQCD predictions for the moments. 

In model (|47p large momentum tail which can not be calculated reliably in the framework of potential model is 
removed. Due to this the contribution of relativistic motion of quark-antiquark pair inside the rjc meson to LCWF 
becomes strongly suppressed. As we have already noted from the point of conformal expansion ^ this suppression 
can be achieved by fine tuning of the coefficients a„. But if this fine tuning takes place at some scale ( in paper 
this happens at scale /i ~ /ig) it will be broken at larger scales where effective LCWF (|47p is applied. 

To illustrate this point let us consider the calculation of the cross section of the process e+e" J/vf -I- rjc at energy 
^/s = 10.6 GeV with the LCWF proposed in our paper. Unfortunately to calculate the amplitude of the process 
e^e^ J 1^1 + lie one must know not only leading twist wave functions but also next-to-leading twist wave functions 
of rjc and J /'^ mesons (see papers our calculation the following model for LCWF of r]c and J/^* mesons 

will be used: 

(50) 

where 0^''(^) is the asymptotic LCWF that corresponds to the function (f)i{£,,fJ.), the function (j>{£,,n) was defined by 
equation ([38|) . One can expect that at scale /i ^ /xq this model is rather good for LCWFs for r/c and J/^ mesons 



^ Note, however, that the NRQCD-formahsm order-da corrections are different from the hght cone formalism order-«s corrections. 
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since at this scale rjc and J/^ mesons have equal moments up to the corrections of higher order in relative velocity 
expansion. At larger scales the model is not so good since the evolutions of LCWFs n) with /i are different. But 
we believe that it is sufficient to apply this model to the estimation of the cross section of the process e+e^ J/'i+rjc- 
It should be noted here that in model (150)1 large relative momentum tail of the LCWFs at scale /i ~ /j,o is strongly 
suppressed as it was required in paper 22]. The calculation of the cross section of the process e~^e~ — *■ J/"^ + rjc [3 
with the wave function taken at scale /i /io gives the result 

cr(e+e- ^ J/* + 77c) 12.1 fb. (51) 

This result is larger than cr = 8.19 fb obtained in paper since LCWF ([55]) is wider than l[Tf|) . 

If the same calculation is done with LCWF taken at scale approximately equal to the momentum running through 
gluon propagator (see Fig. 3) \/k^ — \/s/2 = 5 GeV than we have 

cr(e+e^ -> J/^- + 77c) = 25.1 fb. (52) 



LCWF proposed by Ma and Si(MS) in paper (3C 

Last light cone distribution amplitude that will be considered in our paper is the LCWF proposed in paper [30|. 
The authors of this paper build LCWF as a series in as 

0Bs(C, m) = 0°(e) + <f>\i, f^) + 0{al). (53) 

At lea ding order approximation the function (jp ^ S{^) was taken. Next-to-leading function 4>^{£,,fJ.) can be found in 
paper [3Q|. Due to the singular function 0^(^,^) and infinitely narrow (/)o(0 at scale /x ~ //q the whole LCWF ([53]) 
has negative moments 

(Dms = / Cc^Msit M = "^c) < 0. (54) 

For this reason we do not show the moments of this LCWF in Table II. It should be noted here that as it was shown 
above LCWF of T^c-meson is rather wide and we don't think that it is possible to consider it as an infinitely narrow 
object at leading order approximation. 



VIII. CONCLUSION. 



In this paper we have calculated the moments of leading twist light cone wave function (LCWF) of i]c mesons 
within three approaches. In the first approach we have applied BuchmuUer-Tye and Cornell potential models to 
the calculation of the moments of LCWF. In the second approach the moments of LCWF were calculated in the 
framework of NRQCD. In the third approach the method QCD sum rules was applied to the calculation of the 
moments. The results obtained within these three approaches are in good agreement with each other for the second 
moment (^^). There is a little disagreement between the predictions for the fourth moment {^'^). The disagreement 
between the approaches becomes dramatic for the sixth moment (5^). The cause of this disagreement consists in 
the considerable contribution of relativistic motion of quark- ant iquark pair inside rjc meson to higher moments which 
cannot be regarded reliably in the framework of potential models and NRQCD. The approach based on QCD sum 
rules is more reliable, especially for higher moments since it does not consider ?7c-meson as a nonrelativistic object. 
The main problem of QCD sum rules is that since there is no expressions of radiative corrections to sum rules one 
does not know the size these corrections. However one can expect that QCD radiative corrections will not change our 
results for the moments n = 2 and n = 4 dramatically. As to the sixth moment, the contribution the QCD radiative 
corrections in this case may be important. 

Based on the values of the moments obtained within QCD sum rules we have proposed the model of LCWF. This 
model has some interesting properties: 

1. Due to the evolution ([5]) the accuracy of the moments obtained within the model (|38p improves as the scale rises. 
For instance, if the error in determination of the moment (^^) is 10% at scale ^ = iiq — 1.2 GeV, at scale ^ — 10 GeV 
the error is 4%. For higher moments the improvement of the accuracy is even better and there is no error at all at 
infinite scale /i = oo. The improvement of the accuracy allows us to expect that the model (|38p will be rather good 
even after inclusion of the radiative corrections. 

2. At scale fi ^ the LCWF can be considered as wave function of nonrelativistic object with characteristic width 
~ ~ 0.3. Due to the evolution, at larger scales relativistic tail appears. This tail cannot be considered in the 
framework of NRQCD and at these scales, strictly speaking, rjc meson is not a nonrelativistic object. 
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3. We have found that due to the presence of high momentum tail in the LCWF at scales n > no there is violation 
of velocity scaling rules obtained in paper [l[. More exactly, for any scale fi > there is a number hq for which the 
moments (C"),n > no violate NRQCD velocity scaling rules. 

Actually the last two properties are properties of real LCWF of rjc meson. 

In the last section of our paper we have compared the model of LCWF (|38p with other model proposed in literature. 
Thus we have considered the following models: 

1. Asymptotic LCWF (|43p . Asymptotic LCWF is wider than LCWF ([55)1 . From this comparison we have 
concluded that although rjc meson at scale ^ fiQ can be considered as a nonrelativistic object but the width of this 
meson is rather large and this approximation is not very good. At larger scales LCWF of 77c is wide and r/c cannot be 
considered as a nonrelativistic object. 

2. The LCWF proposed by Bondar and Chernyak in paper |l3l |. The model proposed in paper [13] is in 
good agreement with our model of wave function taken at scale ^ = 10 GeV. 

3. Effective LCWF proposed by Bodwin, Kang and Lee in paper [22,] . Effective LCWF of rjc meson 
proposed in paper is considerably narrower than our model of LCWF. Moreover in calculation of the cross section 
of the process e~^e~ J /'^ + 77c the authors of paper did not regard the evolution of their wave function from 
scale /i ~ ^0 to characteristic scale of this process ~ y/s = 10.6 GeV. We have proposed simple model (|5D|) to estimate 
this effect. Within this model we have shown that the effect is important and good agreement with experimental 
results can be achieved. 

4. LCWF proposed by Ma and Si (MS) in paper (30l |. In this model LCWF is built as a series in ag- At 
first approximation the author took ~ <5(^). Due to the singular next-to- leading order function and infinitely narrow 
leading order function LCWF proposed in paper |33| has negative moments. 
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